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Abstract. 

Invisible plasma content in blazar jets such as protons and/or thermal 
electron-positron (e^) pairs is explored through combined arguments of dynam- 
ical and radiative processes. By comparing physical quantities required by the 
internal shock model with those obtained through the observed radio-to-gamma- 
ray spectra for Mrk 421, we find the existence of a copious amount of invisible 
plasma in the jet. We speculate that the blazar sequence could arise from vari- 
ations of total amount and/or blending ratio of pair and electron-proton 
plasma. 



1. Introduction 

The plasma content of relativistic jets is not easily constrained by observations 
since the emission is dominated by that from non-thermal electrons and it is 
difficult to directly constrain thermal plasma. Hence, the plasma composition 
in AGN jets, whether normal proton-electron (e/p) plasma or electron-positron 
pairs (e^) is a dominant composition, is still a matter of debate (e.g., Reynolds 
et al. 1996; Wardle et al. 1998; Hirotani et al. 1999; Sikora & Madejski 
2000). This problem prevents us from estimating the total mass and energy flux 
ejected from a central engine. In order to constrain the invisible matter in the 
jets, dynamical considerations are indispensable. In Kino & Takahara (2004), 
we proposed a new procedure to constrain the invisible plasma component in 
classical FR II radio sources. We used the fact that the mass and energy densities 
of the sum of thermal and non-thermal (NT hereafter) particles are larger than 
those of non-thermal electrons which are determined by observations. We apply 
the same technique to the inner core jets of AGNs (i.e., blazars) based on the 
standard internal shock model (Kino & Takahara 2008). 



2. Mass density of invisible plasma 

The methodology of constraining the invisible plasma content in the emission 
region is as follows. A lower limit to the total mass density (sum of non-thermal 
electrons and invisible plasma) is restricted by the definition that the mass 
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density of total plasma (p) should be smaller than that of the non-thermal 
electrons (p^ T ). Similarly, a lower limit to the total internal energy density 
(e) is limited by the definition that e is smaller than that of the NT elec- 
trons (e^ T ). Using a relativistic shock jump conditions, we can estimate the 
values of p and e (see Kino & Takahara 2008). The mass density of non- 
thermal electrons of TeV blazars can be well estimated by the fitting of a 
multi-frequency spectrum with synchrotron self-Compton (SSC) model (Kino 
et al. 2002, hereafter KTK). The number density of non-thermal electrons 
n^ T has been estimated as n^ T = f°° . n e ('j e )d'y e , while e^ T is given by 

e^ T = (7 e )n^ T m e c 2 , where n e (^ e ) and (j e ) are the energy spectrum and the 
average Lorentz factor of NT electrons, respectively. By a detailed comparison 
of the SSC model with the observed broadband spectrum of Mrk 421, we ob- 
tained n^ T as n^ T ~ 11 x (7 e ,min/10)~ ' 6 cm" 3 . Here we adopt the index of 
injected electrons for Mrk 421 as s = 1.6 (e.g., Kirk &: Duffy 1999). As for the 
average energy of NT electrons, we obtained 

e^/nf 1 = (7 e )m e c 2 ~ 3.1 x 10 2 m e c 2 , (1) 

where (7 e ) is the average Lorentz factor of NT electrons. Since for s = 1.6, 
electrons near the cooling break energy 7 ej br ~ 10 carry most part of the kinetic 
energy and e^ T has a weak dependence on the minimum Lorentz factor 7e jiri i n . 
In Fig. 1, we show that the case of 7 e ,min ra 10 4 is ruled out for Mrk 421 since 
the case does not fit the EGRET data (KTK). Therefore, Eq. (pQ) is justified in 
any case for Mrk 421. 

The condition of e^ T /e = ('y e )p^ T /T Te \p < 1 can be rewritten as 

J_ (7e)„, 3.1xl0 2 

r rcl - r rel ' {Z > 

where r re i is the relative bulk Lorentz factor between an upstream and a down- 
stream of the shock. From this, we directly see that the invisible mass density 
at least about 100 times the rest mass density of NT electrons is definitely re- 
quired in the framework of internal shock model in which r re i < afew takes 
place. In other words, a loading of baryons and/or a thermal pair plasma is 
expected. As for the upper limit, we use the constraint that bremsstrahlung 
emission from thermal electron (and positron) component should not exceed the 
observed 7-ray emission. We can thus bracket the amount of total mass density 
in the emission region from below and above. In this way we apply this method 
to the TeV blazar Mrk 421. 

In Fig. 2 we show the resultant p/p^ T plotted versus r r /T s in the cases of 
"equal p" and "equal mass" (m = pTA) where r r , r s , and A denote the Lorentz 
factors of rapid and slow shell, and the thickness of the shell measured in ISM 
frame, respectively. We impose the condition A r /A s = 1 (e.g., Spada et al. 
2001). As for the plasma content, we examined the two simple cases. One is the 
case of the jet with pure pair plasma content, whilst the other is the jet made 
of pure e/p plasma. The qualitative features are the same for these two cases, 
although different in quantitative detail. Central to the results is that a large 
amount of mass density of invisible plasma is clearly required in the emission 
region. As the value of r r /r s increases, the value of r re i becomes larger and the 
lower limit on the invisible mass density (p/p^ T ) reduces. 
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Figure 1. One-zone SSC model spectra for the steady state emission of Mrk 
421 (KTK). The thick solid line in the left one shows the best fit spectrum. 
In the right one, we examine the case of 7 e , m i n = 10 4 (the thick solid line). 
The case of 7 e , m in = 10 4 is ruled out because it cannot fit the spectrum well. 



3. Discussion 



Let us discuss the origin of blazar sequence (Fossati et al. 1998) which is tightly 
connected to the nature of the central engine. In Fossati et al. (1998) they com- 
puted average spectral energy distributions from radio to gamma-rays for com- 
plete sample of blazars. They claimed that the continuous sequence of properties 
may be controlled by a single parameter, related to the bolometric luminosity. 
Here we enlighten another new ingredient for the origin of sequence. Flat spec- 
trum radio quasars (FSRQs) have the order of magnitude larger kinetic power 
and larger energy densitity of the surrounding external radiation field than BL 
Lacs (e.g., Sikora et al. 1997). Hence the leptonic components in FSRQs ejecta 
undergo severe radiation drag in dense external radiation fields (Sikora & Wilson 
1981; Phinney 1982). However, the bulk Lorentz factors in FSRQs are slightly 
larger than the ones in TeV blazars (Jorstad et al. 2001; Kellermann et al. 2004) 
in spite of being subject to much stronger radiation drag (e.g., Kubo et al. 1998; 
Spada et al. 2001). In order to realize the larger kinetic powers and larger bulk 
Lorentz factors against the strong radiation drag, we may take a new conjecture 
that a larger baryon loading may occur for FSRQs. Summing up, not only the 
strength of the external radiation field but also the total amount and/or blend- 
ing ratio of pair and e/p could be new key quantities to explore the origin of 
the continuous blazar sequence. 
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Figure 2. The allowed regions (the gray shaded regions) of p/p^ T : the 
amount of mass density of total plasma normalized by that of non-thermal 
electrons for "equal p" (left) and "equal m" (right). Horizontal axis shows 
the ratio of the Lorentz factor of a rapid shell to a slow one. 



References 

Fossati, G. Maraschi, L. Cclotti, A. Comastri, A. & Ghisellini, G. 1998, MNRAS, 299, 
433 

Hirotani, K. Iguchi, S. Kimura, M. & Wajima, K. 1999, PASJ, 51, 263 

Jorstad, S. G. Marscher, A. P. Mattox, J. R. Wehrle, A. E. Bloom, S. D. & Yurchenko, 

A. V. 2001, ApJS, 134, 181 
Kellermann, K. I. ct al., 2004, ApJ, 609, 539 
Kino, M. Takahara, F. & Kusunose, M. 2002, ApJ, 564, 97 (KTK) 
Kino, M. & Takahara, F. 2004, MNRAS, 349, 336 
Kino, M. & Takahara, F. 2008, MNRAS, 383, 713 

Kirk, J. G. & Duffy, P. 1999, Journal of Physics G Nuclear Physics, 25, 163 

Kubo, H. Takahashi, T. Madejski, G. Tashiro, M. Makino, F. Inouc, S. & Takahara, F. 

1998, ApJ, 504, 693 
Phinney, E. S. 1982, MNRAS, 198, 1109 

Reynolds, C. S. Fabian, A. C. Celotti, A. & Rees, M. J. 1996, MNRAS, 283, 873 
Sikora, M. & Wilson, D. B. 1981, MNRAS, 197, 529 

Sikora, M. Madejski, G. Moderski, R. & Poutanen, J. 1997, ApJ, 484, 108 
Sikora, M. & Madejski, G. 2000, ApJ, 534, 109 

Spada, M. Ghisellini, G. Lazzati, D. & Celotti, A. 2001, MNRAS, 325, 1559 
Wardle, J. F. C. Homan, D. C. Ojha, R. & Roberts, D. H. 1998, Nature, 395, 457 



